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ABSTRACT 

New integral field spectroscopy has been obtained for IZwI8, the nearby lowest-metallicity galaxy 
considered our best local analog of systems forming at high-z. Here we report the spatially resolved 
spectral map of the nebular HeHA4686 emission in IZwl8, from which we derived for the first time 
its total Hell-ionizing flux. Nebular Hell emission implies the existence of a hard radiation field. 
Hell-emitters are observed to be more frequent among high-z galaxies than for local objects. So 
investigating the Hell-ionizing source(s) in IZwI8 may reveal the ionization processes at high-z. Hell 
emission in star-forming galaxies, has been suggested to be mainly associated with Wolf-Rayet stars 
(WRs), but WRs cannot satisfactorily explain the Hell-ionization at all times, in particular at lowest 
metallicities. Shocks from supernova remnants, or X-ray binaries, have been proposed as additional 
potential sources of Hell-ionizing photons. Our data indicate that conventional Hell-ionizing sources 
(WRs, shocks. X-ray binaries) are not sufficient to explain the observed nebular HeHA4686 emission in 
IZwI8. We find that the Hell-ionizing radiation expected from models for either low-metallicity super- 
massive O stars or rotating metal-free stars could account for the Hell-ionization budget measured, 
while only the latter models could explain the highest values of HeHA4686/H/3 observed. The presence 
of such peculiar stars in IZwl8 is suggestive and further investigation in this regard is needed. This 
letter highlights that some of the clues of the early Universe can be found here in our cosmic backyard. 

Subject headings: galaxies: dwarf — galaxies: individual (IZwl8) — galaxies: ISM — galaxies: stellar 
content — ISM: lines and bands 

1. INTRODUCTION 

Hell recombination emission indicates the presence 
of very hard ionizing radiation with photon energies > 

54eV. Star-forming galaxies with lower metallicities tend 
to have larger nebular HeHA4686 line intensities com¬ 
pared to those with higher metallicities (e.g. Guseva et 
al. 2000; Schaerer 2003). While nebular Hell emission 
has been observed in some local low metallicity (Z) star- 
bursts (e.g. Schaerer et al. 1999; Guseva et al. 2000; 

Kehrig et al. 2004; Thuan & Izotov 2005), Hell-emitters 
are apparently more frequent among high-redshift (z) 
galaxies than for local objects. Recent work has found 
that > 3% of the global galaxy population at z~3 show 
narrow Hell lines (Gassata et al. 2013), while this num¬ 
ber is much lower at z^O (Kehrig et al. 2011). The 
Hell lines have been suggested as a good tracer of Pop¬ 
ulation HI stars (PopHI-stars; the first very hot metal- 
free stars) in high-z galaxies (e.g. Schaerer 2003, 2008). 

These stars, which should produce a large amount of hard 
ionizing radiation, are believed to have contributed sig¬ 
nificantly to the Universe’s reionization, a challenging 
subject in contemporary cosmology (e.g. Bromm 2013). 

Before interpreting the emission-line spectra of distant 
star-forming galaxies, it is crucial to understand the for¬ 
mation of high-ionization lines in the nearby Universe. 

The ideal place to perform this study is in extremely 
metal-poor nearby galaxies with nebular Hell emission, 
which are the natural local counterparts of distant Hell- 
emitters. 


In this regard, we have been carrying out a programme 
to investigate nearby low-Z starburst systems using the 
integral field spectroscopy technique (e.g. Kehrig et al. 
2008,2013; Perez-Montero et al. 2011,2013). As a part 
of this programme, we have recently obtained new deep 
integral field spectroscopic (IFS) data of lZwl8. This 
is a nearby (D = 18.2 Mpc; Aloisi et al. 2007)0 Hll 
galaxy, well known for its extremely low Z ^ 1/32 so¬ 
lar (e.g. Vilchez & Iglesias-Paramo 1998), which keeps 
lZwl8 among the three most metal-poor galaxies known 
in the local Universe (e.g. Thuan et al. 2004). Its ob¬ 
servational characteristics make lZwl8 an excellent local 
analogue of primeval systems (see e.g. Lebouteiller et al. 
2013 and references therein). 

The presence of the nebular HeHA4686 line in the spec¬ 
trum of lZwl8 has been reported before, though the pre¬ 
cise location and extension of this HeHA4686 emission 
is not known (e.g. Garnett et al. 1991; Izotov et al. 
1997; Legrand et al. 1997; Vilchez & Iglesias-Paramo 
1998). Our unique IFS data unveil for the first time 
the entire HeHA4686-emitting region and its structure in 
lZwl8 (see section 3.2). 

Despite various attempts to explain it, the origin of 
the nebular Hell emission in Hll galaxies/regions still 
remains difficult to understand in many cases; several 
potential mechanisms (e.g. hot Wolf-Rayet (WR) stars, 
shocks from supernovae remnants, X-ray sources) are 

^ A distance of 18.2 Mpc is assumed in this work 
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Fig. 1.— Colour-composite image of IZwlS (blue = Ha from 
Palomar, green = far-UV/GALEX, red = SDSS r’). The box rep¬ 
resents the FOV (16”xl6”) of the PMAS spectrograph over the 
galaxy main body and the extended Ha halo. The PMAS FOV 
is centered on the coordinates RA(J2000.0)= 09^:34™':02®.2 and 
Dec.(J2000.0)= -|-55°:14':25" 


proposed to account (in part or fully) for the Hell ioniza¬ 
tion in these objects (e.g. Garnett et al. 1991; Schaerer 
1996; Dopita & Sutherland 1996; Cerviho, Mas-Hesse & 
Kunth 2002; Thuan & Izotov 2005; Kehrig et al. 2011; 
Shirazi & Brinchmann 2012). Though hot WRs have 
been suggested before as the source of Hell-ionizing pho¬ 
tons in IZwlS (e.g. Izotov et al. 1997; de Mello et al. 
1998), the main mechanism powering the nebular Hell 
emission in this galaxy is still an open issue. 

In this letter, using new IFS data, we derive for the first 
time the total HeHA4686-ionizing flux in IZwl8 and pro¬ 
vide new clues to constrain the sources of high-ionization. 

2. INTEGRAL FIELD SPECTROSCOPIC DATA 

We carried out new IFS observations of IZwl8 us¬ 
ing the Potsdam Multi-Aperture Spectrophotometer 
(PMAS; Roth et al. 2005) on the 3.5 m telescope at 
the Calar Alto Observatory (Almeria, Spain). The data 
were taken in 2012 December with a typical seeing of 
1”. Each spaxel has a spatial sampling of 1” x 1” on the 
sky resulting in a field-of-view (FOV) of 16” xl6” (~ 1.4 
kpc X 1.4 kpc on IZwl8; see Figure[T]). One pointing of 
IZwlS, encompassing its main body which hosts the two 
brightest stellar clusters (referred to as the NW and SE 
knots), was taken during 2.5 hours integration split into 
six exposures of 1500 seconds each. We used the V500 
grating which covers from ~ 3640 to 7200 A and pro¬ 
vides a linear dispersion of 2 A/pixel, and full width 
at half-maximum effective spectral resolution of ^ 3.6 A. 
Calibration images (exposures of standard star, arc and 
continuum lamps) were also obtained. The data reduc¬ 
tion was performed following the procedure described in 
Kehrig et al. (2013). 

3. RESULTS 

3.1. Emission-line flux maps 


The emission-line fluxes were measured using the 
IRAJ0 task SPLOT. The flux of each emission line was 
derived by integrating between two points given by the 
position of a local continuum placed by eye. For each 
line, this procedure was repeated several times by vary¬ 
ing the local continuum position. The final flux of each 
line and its associated uncertainty were assumed to be 
the average and standard deviation of the independent, 
repeated measurements (e.g. Kehrig et al. 2006). 

We used our own IDL scripts to create the emission¬ 
line maps presented in Figure [21 The spaxels where we 
measure HeHA4686 are indicated with pluses. 

3.2. The spatially resolved HeIIX4686-emitting region 

An extended HeHA4686-emitting region with a diam¬ 
eter of w 5” (« 440 pc) is revealed from our IFS data 
(see Figured]). The narrow line profile for the HeHA4686 
emission and its spatial extent are evidence for its nebu¬ 
lar nature. The spectra of certain Of stars may exhibit 
somewhat narrow HeHA4686 lines (e.g. Massey et al. 
2004), so Hell emission observed in starburst galaxies 
could be thought to arise in the atmospheres of such stars 
(e.g. Bergeron 1977). If these Of stars were present in 
appreciable numbers in IZwl8, broad Ha emission would 
however be expected (Massey et al. 2004). Thus, the 
non-detection of this feature in our spectra supports the 
nebular origin for the HeHA4686 line in IZwl8. 

By adding the emission from the spaxels showing 
HeIIA4686 (see Figured]), we created the one-dimensional 
spectrum for the HeHA4686 region (see Figured]). Us¬ 
ing this spectrum, we obtained a very small logarith¬ 
mic extinction coefficient c(H/3) = 0.08 ± 0.02 from 
the observed ratio Ha/H,5 = 2.92 ± 0.04, assuming an 
intrinsic case B recombination Ha/H/3 = 2.75 (Oster- 
brock & Ferland 2006, OF06). Using PYNEB (Lurid- 
iana, Morisset & Shaw 2015), we derived the electron 
temperature, Tg = (2.18 ± 0.05) x 10^ K, and density, 
Ue < 100 cm-3, from the [OHI]A4363/[OHI]AA4959,5007 
and [SH]A6717/[SH]A6731 ratios, respectively. The inte¬ 
grated flux of the HeHA4686 line, corrected for redden¬ 
ing, is (2.84 ± 0.18) X10“^® erg s“^ cm“^ which trans¬ 
lates to a HeHA4686 luminosity of Tneiixiese = (1.12 
± 0.07)X10^® erg s“^. The corresponding Hell ioniz¬ 
ing photon flux, Q(HeH)obs = (1.33 ± 0.08) x 10®° pho¬ 
ton s“^, was derived from the measured THeiiXAese using 
the relation Q(HeH) = Liye//A4686/[j(A4686)/aB(HeH)] 
(assuming case B recombination, and Tg = 2x10^ K; 
OF06). We checked that different CLOUDY models 
(Ferland et al. 2013) computed for the typical Ue and 
Te in the HeHA4686 region, and considering several ef¬ 
fective temperatures and geometries with no dust, give 
a ratio L(HeHA4686)/Q(HeH) which agrees with the as¬ 
sumed ratio (OF06) within 10%. The total Q(HeH)o6s, 
a quantity not reported before for IZwl8, will allow us 
to constrain possible ionizing sources of Hell in IZwlS. 

4. DISCUSSION 
4.1. Ionizing sources of Hell 

^ IRAF is distributed by the National Optical Astronomical Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 2. — Emission-line flux maps of IZwl8. Maps are displayed in logarithmic scale and the fluxes are in units of erg s~^ cm“^; the area 
of each spaxel is 1 arcsec^ on the sky. Top row. Ha and HeIIA4686 maps. Middle row: for display purposes, the maps of Ha and HeHA4686 
are presented as color filled contour plots and were smoothed using bilinear interpolation. Isocontours of the Ha emission line flux are 
shown overplotted as reference. Bottom row: [OIJA6300 and [SH]AA6717-I-6731 maps. The spaxels where we detect nebular HeHA4686 
are marked with pluses on the maps of Ha (top row), [OIJA6300 and [SH]AA6717-I-6731 (bottom row). The spaxels with no measurement 
available are left blank. 
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Fig. 3. — Portion (~ 4425 - 5100 A) of the integrated spectrum 
of the HeIIA4686-emitting region of IZwl8. 

One widely favored mechanism for Hell ionization in 
HII galaxies involves hot WRs (e.g. Schaerer 1996). 
Nevertheless, it has been demonstrated that nebular 
HeIIA4686 does not appear to be always associated 
with WRs, as it is the case of the Hell nebulae LMC- 
N44C, LMC-N159F and M33-BCLMP651, among others 
(Kehrig et al. 2011 and references therein). This in¬ 
dicates that WRs cannot explain Hell ionization in all 
cases, particularly at low-Z (e.g. Guseva et al. 2000; 
Kehrig et al. 2008,2013; Schaerer 2003; Shirazi & Brinch- 
mann 2012). Besides, it has been shown that, for a sam¬ 
ple of HeHA4686-emitting star-forming galaxies, current 
models of massive stars predict HeHA4686 emission, and 
HeHA4686/H/3 ratios, only for Z > 0.20 Zq instantaneous 
bursts (Shirazi & Brinchmann 2012). In the particu¬ 
lar case of IZwl 8 , previous work has claimed that the 
ratio HeHA4686/H/3 could be reproduced using highly 
density-bounded photoionization models while underpre¬ 
dicting the electron temperature measurement (Stasihska 
& Schaerer 1999); these models have been challenged by 
Pequignot (2008). 

Faint broad emission signatures, attributable to WRs, 
are observed in the spectrum of IZwl 8 despite its low-Z 
(e.g. Izotov et al. 1997, Legrand et al. 1997, Brown et al. 
2002). A comprehensive study of WRs in IZwl 8 , using 
UV STIS spectroscopy, revealed signatures of carbon- 
type WRs (WC) in two clusters: one in the NW star¬ 
forming region and a second one on the outskirts of this 
region (Brown et al. 2002). Here we deal with the NW 
one, since it is in the HeHA4686-emitting region (see Fig¬ 
ure [2]). The CIVA1550 flux measured from the NW WR 
cluster (Brown et al. 2002) gives a CIVA1550 luminosity 
of Li 55 o = 4.67x10^^ erg s“^. Taking the L 1550 luminos¬ 
ity of the metal-poor early-type WC (WCE) model by 
Crowther & Hadfield (2006, CH06), which mimics a sin¬ 
gle WCE star in IZwl 8 , this implies ^ 9 IZwl 8 -like WCE 
stars present in the NW cluster. From these 9 WCE stars 
a total flux Q(HeH) =2.8 xlO^® photon s“^ is expected 
[assuming Q(HeH) = 10^^ ® photon s“^ for one IZwl 8 - 
like WCE; CH06], i.e. about 48 times lower than the 
Q(HeH)of,s=(1.33±0.08)xl0®° photon s“^, derived from 
our data (see section 3.2). 

Based on the Hell-ionizing flux expected from these 
IZwl 8 -like WRs, a very large WR population is re¬ 
quired to explain the Hell-ionization budget measured; 
for instance, taking the Q(HeH)= lO"^^® photon s“^ 
for one IZwi 8 -like WCE (CH06), the number of these 
WCEs needed to explain our derived Q(HeH)obs = 
(1.33±0.08 )x 10'5° photon s'^ would be > 400. In prin¬ 
ciple, the presence of hundreds of WRs in IZwl 8 should 


not be discarded on the basis of empirical arguments for 
reduced WR line luminosity at low-Z (CH06). However, 
assuming a Salpeter initial mass function (IMF; Salpeter 
1955; M„p=15OM0) and the initial mass needed for a star 
to certainly become a WC (Maeder & Maynet 2005), a 
cluster with > 8 times the total stellar mass of the NW 
region (M*_vrvv = 2.9x 10® Mq from Stasihska & Schaerer 
(1999) scaled to 18.2 Mpc distance) is required to pro¬ 
vide > 400 WC in IZwI8. Also, such a high number 
of WRs is not supported by state-of-the-art stellar evo¬ 
lutionary models for single (rotating and non-rotating) 
massive stars in metal-poor environments (Leitherer et 
al. 2014). Furthermore, given the decrease in the ra¬ 
tio of WR/0 stars with decreasing metallicity, shown by 
observations and theoretical models (Maeder & Meynet 
2012 and references therein), such a large number of WRs 
appears clearly unreasonable considering the extremely 
low-Z and the O star content of IZwl8 (CH06). All this 
suggests that WRs are not the sole responsible of the 
HeHA4686 emission in IZwl8. 

The binary channel in massive star evolution is sug¬ 
gested to increase the WR population (e.g. Eldridge et 
al. 2008), but the WR population in Local Group galax¬ 
ies does not show an increased binary rate at lower-Z 
(e.g. Foellmi et al. 2003; Neugent & Massey 2014); thus 
the binary channel does not seem to favour the forma¬ 
tion of WRs at lower-Z, in contrast to what we need. 
Nevertheless, we should bear in mind the possible uncer¬ 
tainties yet unsolved in the models (Maeder & Meynet 
2012). Further investigation awaits the calculation of 
evolutionary models for binary stars at very low-Z. 

As mentioned before, nebular HeHA4686 emission ob¬ 
served in the spectra of HII galaxies and extragalactic 
HII regions has often been attributed also to shocks and 
X-ray sources (e.g. Pakull & Mirioni 2002; Garnett et 
al. 1991; Thuan & Izotov 2005). In the following we dis¬ 
cuss on these two candidate sources for Hell ionization 
in IZwl8. 

The X-ray emission from IZwl8 is dominated by a sin¬ 
gle X-ray binary apparently located in the field of the 
NW knot (Thuan et al. 2004). We have computed a 
CLOUDY photoionization model (Ferland et al. 2013) 
using as input a power-law spectral energy distribution 
(SED) with the same X-ray luminosity, column density 
and slope that have been reported for IZwl8 (Thuan et 
al. 2004). This CLOUDY model gives as an output a 
HeHA4686 luminosity ljHeii \4686 = 10®®'^ erg s“^ which 
is ^ 100 times lower than the L//e//A 4686 measured. This 
result rules out the X-ray binary as the main source of 
Hell ionizing photons in IZwl8. We note here that the 
emission from X-ray ionized nebulae has been success¬ 
fully reproduced by CLOUDY models before (e.g. Pakull 
& Mirioni 2002). 

Guided by the existence of HeHA4686 emission as¬ 
sociated with supernova remnants (e.g. Kehrig et al. 
2011) we explored the conjecture that the HeHA4686 
region represents such a shock-ionized nebula. The 
[OIJA6300 line, often strong in remnants, has been fre¬ 
quently used as a sensitive shock-emission test (e.g. Skill- 
man 1985). We checked that, in fact, most of the 
[OIJA6300 emission in IZwl8 is concentrated on the SE 
knot (see the [OI]A6300 map in Eigure[2]), with only 12% 
of the HeHA4686-emitting spectra showing [OI]A6300 
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flux above the 3cr detection limit. Additionally, we find 
no evidence for [SII] enhancement (a usual sign of shock 
excitation; e.g. Dopita & Sutherland 1996) associated to 
the location of the HeIIA4686 region (see the [SII] map 
in Figure 12). Therefore, the HeIIA4686-emitting zone in 
IZwl 8 is unlikely to be produced by shocks. 

4.2. Peculiar very hot stars in IZwl8 ? 

Our new observations have allowed us to empirically 
demonstrate why conventional Hell-ionizing sources (e.g. 
WRs, shocks. X-ray binaries) cannot account for the to¬ 
tal Hell-ionization budget in IZwl 8 . What could the 
nebular HeHA4686 emission in IZwl 8 originate from ? 

We have also explored the possibilty of very massive, 
metal-poor O stars to account for our observations of 
IZwl 8 . Using current wind models of very massive O 
stars at low-Z we can derive their Hell-ionizing fluxes 
(Kudritzki 2002, K02). According to the hottest mod¬ 
els (re=60,000 K), between 10-20 super-massive stars 
with 300Mq [with Q(HeH) Ri (0.70-1.4) x 10^® photon s“^ 
each] would be enough to explain the derived Q (Hell) obs 
budget. Very massive stars of up to 3 OOM 0 were claimed 
to exist in the LMC-R136 cluster (Crowther et al. 2010); 
however the existence of such super-massive 300Mq stars 
remains heavily debated (Vink 2014). Additionally, as¬ 
suming a Salpeter IMF, 10-20 stars with 290 < M^/Mq 
< 310 would imply a cluster mass of 10-20 x 
We should bear in mind that an extrapolation of the IMF 
predicting 3 OOM 0 stars remains unchecked up to now. If 
we consider instead the 15OM0 star hottest models [with 
Q(HeH) < 1.9x10“^^ photons”^ each, for Z < 1/32 Z 0 ; 
K02], the number of these stars required to explain the 
Q(HeH)ot,s would be > 650. For a Salpeter IMF, 650 
stars with 145 < M*/M 0 < 155 would require a cluster 
mass ^ 200 x Besides the Q(HeH)obs budget, in 

the HeHA4686 region we have measured HeHA4686/H/3 
ratios as high as 0.08. These values appear too big to be 
explained even by the models for the hottest, most metal- 
poor super-massive 3 OOM 0 stars (K02) under ionization- 
bounded conditions. Further constraints to the observa¬ 
tions should await the calculation of new evolutionary 
tracks and SEDs for single O stars at the metallicity of 
IZwl8 including rotation. 

Searches for very metal-poor starbursts and PopHI- 
hosting galaxies have been carried out in the distant Uni¬ 
verse using the Hell lines (e.g. Schaerer 2008; Cassata 
et al. 2013). This search is based on the high effec¬ 
tive temperature for PopHI-stars which will emit a large 
number of photons with energy above 54eV, and also 
on the expected increase of the Hell recombination lines 
with decreasing Z (e.g. Guseva et al. 2000; Schaerer 
2003,2008). Predictions for burst models of different 
metallicities show how their corresponding Q(HeH) can 
increase by up to ^ 10^ when going from Z=10“^ to 
0 (Schaerer 2003). These models cannot explain the 
Q(HeH)of,s when Z > 10“® [for a Salpeter IMF, M„p = 
IOOM 0 ; see table 3 in Schaerer 2003] even assuming that 
the total M*^ 7 vw would come from Hell-ionizing stars. 
So another more speculative possibility, to explain the 
derived Q(HeH)obs could be based on nearly metal-free 
ionizing stars. These stars should ionize Hell via their 
strong UV radiation expected at nearly zero metallicity 
(e.g. Tumlinson & Shull 2000; Schaerer 2003). 


As an approximation of nearly metal-free single stars, 
we have compared our observations with the Hell- 
ionizing radiation expected from state-of-the-art models 
for rotating Z=0 stars (Yoon, Dierks & Langer 2012). 
According to these models, we found that a handful of 
such stars could explain our derived Q(HeH)of,s [e.g. ~ 
8-10 stars with mass Mi„j=15OM0 or ~ 13-15 stars with 
M™=1OOM0; with Q(HeH) « l.dxlO^^ (0.9x10^®) pho¬ 
ton s“^ for each star with Mini =15OM0(1OOM0)]. Ad¬ 
ditionally, we note that the ionizing spectra produced by 
these star models are harder than the ones expected from 
the hottest models of super-massive 3 OOM 0 stars (K02), 
so they would also explain the highest HeHA4686/H/3 
values observed, providing that ionization-bounded con¬ 
ditions are met. While gas in IZwl8 is very metal- 
deficient but not primordial, Lebouteiller et al. (2013) 
have pointed out that the HI envelope of IZwl8 near 
the NW knot contains essentially metal-free gas pock¬ 
ets. These gas pockets could provide the raw mate¬ 
rial for making such nearly metal-free stars. Clearly, in 
this hypothetical scenario, these extremely metal-poor 
stars cannot belong to the NW cluster, which hosts more 
chemically evolved stars. 

5. SUMMARY AND CONCLUDING REMARKS 

This letter reports on new optical IFS observations of 
the nearby dwarf galaxy IZwl8. This is an extremely 
metal-poor system, which is our best local laboratory for 
probing the conditions dominating in distant low-Z star- 
bursts. Our IFS data reveal for the first time the total 
spatial extent (« 440 pc diameter) of the HeHA4686- 
emitting region and corresponding total Hell-ionizing 
photon flux in IZwl8. The metal-poor sensitivity of the 
Hell line is a primary motivation to develop diagnostics 
for unevolved starbursts, and strong nebular Hell emis¬ 
sion is expected to be one of the best signatures of mas¬ 
sive PopHI-stars (e.g. Schaerer 2003,2008). Hell emis¬ 
sion has been observed to be more frequent at higher-z 
than locally (Kehrig et al. 2011; Cassata et al. 2013). 
Thus the analysis of the origin of the HeHA4686 nebu¬ 
lar line in relatively close ionized regions, which can be 
studied in more detail, can yield insight into the ionizing 
sources in the distant Universe. 

Our observations combined with stellar models predic¬ 
tions point out that conventional excitation sources (e.g. 
WRs, shocks, X-ray binaries) cannot convincingly ex¬ 
plain the Hell-ionizing energy budget derived for IZwl8. 
Other mechanisms are probably also at work. If the Hell- 
ionization in IZwl8 is due to stellar sources, these might 
be peculiar very hot stars (perhaps uncommon in local 
starbursts but somewhat more frequent in distant galax¬ 
ies): according to theoretical stellar models, either super- 
massive O stars or nearly metal-free ionizing stars could 
in principle account for the total Q(HeH)of,s of IZwl8. 
However the super-massive O stars scenario would imply 
a cluster mass much higher than the mass of the NW knot 
derived from observations. On the other hand, though 
metal-free gas pockets were previously reported in IZwl8 
(Lebouteiller et al. 2013), we highlight that the existence 
of nearly metal-free ionizing stars is not yet confirmed ob- 
servationally. The work presented here can help in the 
preparation of prospective searches for primeval objects, 
one of the main science drivers for next-generation tele¬ 
scopes (e.g. Bromm 2013). 
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